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SUMMARY 


Hybrids of s, p and d orbitals of the atomic L and M shells have been used in order to describe 
the distribution of the valence electrons of the boron atoms in metal borides of the type MBg. 
Those linear combinations of the hybrid orbitals, which form a basis for the different irreducible 
representations of the symmetry group Op, have been found. Using these orbitals, Bloch wave 
functions have been constructed. In a one-electron approximation the corresponding secular 
equations have been solved numerically by use of the electronic computer Besk. The results show 
that the M shell orbitals are very important for the explanation of the high stability of the erystal. 
The experimentally found conductivities of the crystals of type MB, have been explained theoreti- 
eally on account of conduction bands. 

A preliminary communication of the present investigation has been given at the conference for 
physicists in Uppsala in June 1957 and at the summer school in quantum chemistry at Vala- 
dalen in August 1958. 
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Introduction 


The stable crystal structures of the metal borides of the type MB, have been experi- 
mentally investigated by Kiessling [1]. The notation M in MB, is written for the 
metal atoms, which are either bivalent (Ca, Sr and Ba) or three-valent (Y, La and 
some others of the rare earths). The metal atoms form a cubic framework in the 
crystal. In every cubic cell the boron atoms are grouped in regular octahedrons. The 
unit cell of the crystal is given in Fig. 1 in section 2. 

In an earlier investigation by the present author [2], the MB, crystals were theoreti- 
cally treated in a very simple way by extremely rough valence bond methods. On ac- 
count of the symmetry of the crystal, atomic d orbitals of boron were thought to be 
of some importance. The author also compared the borides MB, to each other for 
different metal atoms. In this connection there was used a rough method of maximum 
charge density in the boron-boron bond directions in comparison with other direc- 
tions. The distance from a boron nucleus to the point of maximum value for the elec- 
tron charge density in these boron—boron bond directions was determined for different 
values of the orbital exponents of the d orbitals. By comparing these distances to the 
experimentally determined distances between the boron atoms for different metal 
atoms in the crystals, it was possible to determine the orbital exponents of the d 
orbitals for the valence electrons of boron in the different metal borides. The result 
showed that the metal borides, ordered in increasing sequence of the ionization poten- 
tials of the metals, also were ordered in increasing sequence with respect to the orbital 
exponents of the d orbitals and also with respect to the stability of the borides in 
decreasing sequence. Thus, the stabilities of the different borides were explained as 
caused by ionization of the metal atoms. 

The MB, crystal has been theoretically investigated also by Longuet-Higgins and 
de V. Roberts [3] in a more accurate way than there has been done in [2], but the 
approximations used in [3] are still rather rough. These authors found, that the boron 
valence electrons together with two electrons of an ionized metal atom in MB, 
form a closed shell structure of 2s and 2p orbitals in each B, octahedron. As is seen 
in section 8, this result is not quite confirmed by the present investigation. 

In the computations of the present work the experimentally determined ionization 
energy and the distances between the atoms of the crystal have been taken for the 
case CaB,, where Ca has been considered to be completely ionized (A = 1 in (38)). The 
corresponding values for the other metal borides of the type MB, are very nearly the 
same. The found results are thus supposed to be rather independent of the special 
kind of metal atoms occurring in the MB, crystal. 

In the present paper an approximate one-electron treatment of the crystal, yielding 
a secular equation, is given in section 1. Bloch orbitals have been used here, con- 
structed by use of non-orthogonal localized orbitals of all the cells. In section 2 
projection operators of the point symmetry group O, for the configuration of the 
atoms inside every cell are used in order to construct basic cellular orbitals for 
MB,. In section 3 there is shown that these orbitals give a natural way of approxi- 
mately factorizing the secular equation. Section 4 gives an approximate one-electron 
Hamiltonian operator for the crystal. Section 5 gives formulae for construction of 
the elements in the split secular equations. In section 6, transformation formulae 
and numerical values for some of the integrals occurring in the secular equations are 
given. The approximations used in the methods of the present investigation are 
discussed in section 7. Section 8 gives the numerical results for the energy bands and 
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the electron distribution. These results are discussed in section 9. The appendix 
contains a derivation of an approximate expression for the Fock operator expressed 
in the representation of a set of atomic orbitals for the different boron atoms in the 
crystal. 

Since the present investigation is very similar to [3], it might be valuable to 
compare the approximations used in the present work by those made in [3]. We can 
distinguish mainly four different points: 

1. In [8] the binding energy integrals have been put proportional to the corre- 
sponding overlap integrals. In the present work the one-electron two-centre energy 
integrals of all the different kinds occurring have been accurately computed, the 
three-centre integrals have been approximately calculated and the two-electron 
Coulomb interaction integrals have been compensated by use of effective charges. 

2. In [3] the one-centre energy integrals have been put equal to each other for 
all the 2s and 2p orbitals. This approximation together with the approximation in 
point 1 cause no energy expressions to occur in the secular equations of [3], which 
only contain overlap integrals. In the present work the solutions of the secular equa- 
tions give the energy bands directly. 

3. In [3] the tight-binding approximation has been used. This means taking into 
account only the 5 nearest B neighbours. In the present work the 24 nearest B neigh- 
bours and the 8 nearest M neighbours are taken into account. Thus also some integrals 
connecting edge-centred cells are included. Numerically, this means a factor of about 
3 for the 2s and 2p overlap elements in the present approximation compared with 
the corresponding elements in the tight-binding approximation. The more distant 
neighbours are very essential on account of their large number. 

4. In [8] only 2s and 2p orbitals have been taken into account. In the present 
work also 3s, 3p and 3d orbitals are considered. Though the contribution of these 
orbitals cannot be calculated with satisfactory accuracy, the numerical results show, 
that these orbitals are very essential for explaining the properties of the crystal. 

The first point above certainly yields the greatest difference between the two 
compared methods. However, also in the present calculations the approximations 
mentioned under the first point are those which are most unreliable. Hence, the 
numerical results cannot be supposed to be quite accurate. In spite of this, the quali- 
tative results found by the present calculations seem to be fairly reliable. 


1. A one-electron approximation in the solid state theory 


Let us consider the Schrédinger wave equation of an arbitrary crystal without any 
external influence. In the Born—Oppenheimer approximation, we can separate the 
electronic wave function from that of the nuclei. Furthermore, we will accept the 
ordinary approximation of considering the nuclei to have fixed positions in the space, 
and we neglect the thermal energy and the zero point vibration energy of the crystal. 
Now, let us assume, for simplicity, that there is an even number, 2N, of electrons in 
the crystal, and suppose that the crystal is completely diamagnetic. Thus, in neglect- 
ing spin-orbit coupling, the total electronic wave function must be an eigenfunction 
of both S? and S, with the eigenvalues zero, where S is the total electronic spin 
operator. Since there are (2N)!/(N!(N + 1)!) different independent spin functions 
of this kind, the total electronic wave function generally is a linear combination of 
these spin functions, where each of them is multiplied by a space wave function of 
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all the electrons. However, in the above-mentioned approximation, the Hamiltonian | 
operator does not affect the spin wave functions. Thus, if we suppose all the electronie¢ | 
states to be doubly occupied, it is sufficient to consider only the space part of the total | 
electronic wave function. Furthermore, we will neglect the correlation between the | 
different electrons. Thus, we can deal with space orbitals of one-electron type only. | 

Now, let us divide the crystal into equivalent parallelopiped unit cells. Following | 
[3], we will use the notation “‘site’’ for the configuration of atoms in the unit cell | 
in what follows. The site can thus be considered as a molecule, which usually does | 
not exist in a free state, but which is stable when forming part of a crystal lattice. 
We introduce the “‘site orbitals’, 


ga(t) =$(r —N), (1). 


which are one-electron wave functions localized to the different sites of the crystal. 
The position vectors of the centres of the sites have been denoted by n. The origin 
of the coordinate system is put at the centre of one of the sites, which we denote 
“the central site’ (n =0). The position of the electron in this coordinate system is 
denoted by r. The site orbitals are generally not orthogonal to each other for different 
n, and they need not necessarily be normalized. According to Léwdin [4], it is possible 
to orthogonalize the site orbitals by use of an expansion of the inverse of the square 
root of the overlap matrix. The elements of this matrix, 


<n[n +m) ={ $F darm dt = Am, (2) 


depend upon the relative position vector, m, of the participating sites. However, 
orthogonalisation is unnecessary here, when we want to construct crystal orbitals, 
yx, of the Bloch type, since these functions are automatically orthogonal to each other. 
The connection between the site orbitals and the crystal orbitals is 


p= (V2 Oars ee a > es" bu (3) 
or daa Vite c= ete ee ee (4) 


Here, the sums over m and n run over all the sites of the crystal, which we consider 
as infinite in all directions. The crystal orbitals are normalized over the so called 
“micro-crystal”, which is the domain for which they are periodic. The micro-crystal 
consists of a three-dimensional infinity of sites, and V = L,L,L, is its infinite volume, 
measured in units of the volume of a single cell. Here, L, are the side lengths of the 
parallelopiped micro-crystal measured in corresponding units. The components of 
k in the space of the reciprocal lattice are 2ah,/L, (v =1,2,3), where h, runs over 
the infinite number of integers of the first zone, ie. —L,/2<h,< L,/2. According 
to the reduced zone scheme [5], we get all the crystal orbitals in this way. In the 
special case, when the site orbitals (4) are orthogonal to each other, (3) turns into 
the ordinary Bloch wave functions [6]. Also when overlap is included, the crystal or- 
bitals (3) have the properties of Bloch functions, i.e. they might be written as exp 
(tkr) times a function, which is periodic over the sites in the crystal. Thus, we have 
the well-known result, that the crystal orbitals (3) form a basis for the one-dimensional 
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irreducible representations of the Abelian group of those translations, which trans- 
form the infinite crystal into itself [6]. 


Neglecting configurational interaction, we now construct an anti-symmetric wave 
function of all the electrons in the crystal as a single Slater determinant of the crystal 


orbitals combined with spin functions. The Schrédinger wave equation then leads 
to the ordinary Hartree-Fock equations, 


F pu= ex Yu; (5) 


where F is the Fock operator and e¢ are one-electron energy eigenvalues. In terms 
of matrix elements in the space spanned by (3), we obtain from (5) 


f= <k| | k> +> {2 (kk | AP | kk — A, (kk || kk}. (6) 


Here A, and A, are certain numerical coefficients, which have been introduced, since 
W. are supposed not to be spin orbitals, but doubly occupied space orbitals. H* and 
H” are the one- and two-electron operators respectively occurring in the total Hamil- 
tonian operator, 


1 
H=H,+>H'+ ay He (7) 
i i 


Introducing the site orbitals (4), it is readily seen that 


wig ee H 
Ck|E'|k> =f ye Di pdt= Se (8) 
Here we have introduced the notation 
Hi=<(m|H|m+n>= [$n A" bara dt (9) 


for the one-electron energy matrix elements in the space spanned by (4). The corre- 
sponding overlap matrix elements, A,, are given by (2). A form quite analogous to 
(8) can be found for the total expression (6) of & (see appendix): 


y e= F, 


&= Syma: 
er he 


(10) 


Here, F, are given by (72) and form matrix elements of the Fock operator F (cf. 
74)). 

in the case, when there are many atoms in a single site, it is sometimes convenient 
to form every site orbital of the type (4) as a linear combination of given atomic 
orbitals, ya, of the same site, n. The index y then refers to a certain atom and a certain 
atomic orbital localized to this atom. This corresponds closely to the ordinary way 
of constructing molecular orbitals. Thus, 


gu =>, Or Pas (11) 
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where ©, are unknown constant coefficients. In introducing (11) into (10), we obtain 
the following expression for the energy parameters: 


N? Ge Cy» Vie 
a ype 
k —— k yt py “| 
3 OnE OL SNA 
ppt 


(12) 


Here we have put 
Kee => eg, Bed (13) 


In (13) there occur atomic orbital integrals containing the Fock energy operator F 
(cf. (75)). The notations Ay”” in (12) stand for expressions in atomic overlap inte- 
grals, which are completely analogous to (13). In the appendix the full expression 
of (13) is given in matrix elements of the one- and two-electron Hamiltonian operators. 

We now suppose that the Fock operator occurring in (5) is replaced by an approxi- 
mate one-electron Hamiltonian operator h, found by assuming a fixed crystal poten- 
tial, which we shall discuss in section 4. Thus, as is well-known, the eigenvalue equa- 
tion (5) is equivalent to the variation problem, 


O& =0, (14) 


where ¢, are one-electron energy eigenvalues closely corresponding to (12). Here 
we choose the coefficients, C,, as variables. Since the ortho-normalisation condition 
automatically is taken into account in (12), there is not needed any subsidiary condi- 
tion to the variation problem (14). 

As is seen in the appendix, Fy” depend upon the coefficients C, and are thus of 
course actually not independent of the variation. However, since it is far from 
practicable to solve the accurate equations, we postulate that hy”” are independent 
of the variation, which corresponds to using a fixed crystal potential. The justifica- 
tion of this approximate attempt depends very markedly on the way of choosing the 
effective energy operator h (cf. section 4). 

On account of this statement, it is now readily verified that (12) and (14) give a 
linear homogeneous system of equations for the coefficients C,. The corresponding 
secular equation, 


Det {hy” —e, Ax” }=0, (15) 


has the roots ex (4 = 1,2,3, ...) for which the linear system of equations can be solved. 
To get the energy bands, the secular equation must be solved for different values of k. 

In the case of many atoms inside every cell, the degree of the secular equation is 
often very large. Then it is advisable to try to factorize this equation. This procedure 
may be facilitated by use of group theory, if the system of atoms of a single site has 
any degree of symmetry. In the following section we shall make an investigation of 
this kind for the case of MB,. 


2. Construction of basic site orbitals for MB, 


; We shall now attack the problem of simplifying the secular equation (15) by apply- 
ing a method given by Wigner [7]. Let us first consider a configuration of atoms, 
which belongs to the symmetry group G. Let R be an element of @, and denote by 
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$7‘ a set of basic orbitals! of the jth irreducible representation. Here j andi can take 


“the values j =1,2,3, ... J and 1=1,2,3,...1,. The matrix elements of R in this 
representation, 


Ru={o™ Rg" dt, 


are coefficients in the expansion 


Rgt= 2 / Bui o". (16) 
These elements satisfy the well-known orthogonality relations, 
SOF pe eS 
es Ri? Rie = — Oy Oi Oxck's 3 (17) 
R L; 


where G@ is the number of elements of the group. On account of (16) and (17), it is 
found that 


+e 
pi SiRER (18) 
R 


is a projection operator, which gives the “‘projection” on ¢’' of any function y which 
can be expanded in the complete function space spanned by ¢’* [7]. 

According to Melvin [8] the projection operator (18) can be factorized in order 
to be more tractable in practical problems. This simplification has, however, not been 
used explicitly in the present investigation, since it is preferable to use sums rather 
than products in the special procedure described below (see Tables 1-5). Factoriza- 
tion has been used only for the “local”? operators, which are referred to the different 
octahedron corners of the central site (see Table 3). 

The constant factor l;/G has no practical importance for the problem of finding 
those basic functions ¢’', which correspond to a certain given arbitrary function. 
Because of that, we will put this factor equal to unity in the following. 

Let us now consider the MB, crystal. Since both the cube and the octahedron are 
transferred into themselves by application of the elements of the octahedral group, 
O,, the whole crystal belongs to this group. The total number of elements in O,, is 
48. In Fig. 1 the axes of rotation for the different elements of the subgroup O of 
O, are given. The notations of the symmetry axes occurring in this figure are the 
same as those used by Melvin [8]. 

Now we interid to find a set of basic site orbitals, 4’(r), of the central site to which 
we refer the group operators, R, in terms of spd hybrid atomic orbitals, localized to 
the six boron atoms of this site. In order to do this, we apply the projection operator 
(18) to a linear combination, ¢, of spd hybrids of the type 


4 = Oy pes + Ce Gone + Cg P2my + C, von. + Cs pas + Cg Gare + Cz P3ry + Cy Pare 
“5 ie P3da =f Cy P3dex Hk Cn P3dzy “ie Ox P3dzy a C13 P3dx2-y2» (19) 


where the (’s are arbitrary constant coefficients and the y's are Slater atomic orbitals 
of the L and M shells, localized to one of the octahedron corners in Fig. 1 (cf. (11)). 
The linear combination goes over all the six corners A, B,... F, thus giving 6 x 13 = 
78 independent coefficients C. In the following we denominate the expansion (19) 


1 In Melvin’s [8] notations these are called symmetry adapted orbitals. 
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Dy/ 


yan 
ee 


Fig. 1. The unit cell in MB,. The different boron atoms of the octahedron corners are denoted 

by A, B,... F. The metal atoms are denoted by M. Lines denoted by D are two-fold symmetry 

axes. Lines denoted by 7’ are three-fold symmetry axes. The coordinate lines are four-fold 

symmetry axes, C,,C, and C,. For the “central site”, O is origin also of the crystal coordinate 
system. 


as the LM approximation. The corresponding expansion finished by the C, term is 
denoted the Z approximation. 

The different elements of the symmetry group O,, transform (19) into “symmetry 
orbitals’? of the same kind as the original one, but localized to the same or to one of 
the other corners of the B, octahedron. When R operates on one of the symmetry 
orbitals the resulting function will be another symmetry orbital. All these symmetry 
orbitals, which we denote by x, can be labelled by four suffixes, 2, bu, vy and o. Here, 
A, « and y are the Eulerian angles of that rotation, which corresponds to R, divided 
by 2/2, and o indicates reflexion. The Eulerian angles are referred to the coordinate 
system of the central site, given in Fig. 1, where the z-axis has been chosen perpen- 
dicular to the line of nodes. Thus, 4 can take the values 0, 1 or 2; jand y can be 0, 1, 2 
or 3. For A=0 and 2, the line of nodes has been chosen in the positive y direction. 


* These are usually called “equivalent orbitals” and must not be confused with the symmetry 
adapted orbitals, ¢’'. 
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Fig. 2. The local coordinate systems for the different corners A, B, ... F of the B, octahedron. 


In this case only two suffixes are needed for the Eulerian angles (see the second 
column in Table 1). The first one of these stands for A and the second one for y + u, 
when A = 0, and »y — w, when A = 2. When R yields no reflexion, o is omitted. In other 
cases o is replaced by a “‘local” operator x. A local operator has the property that 
when it is working on a function localized to the corresponding B atom the result 
will be a function which is also localized to this B atom. The local operator x means 
reflexion in a plane perpendicular to the x-axis of the special coordinate system, which 
belongs to the actual B atom. These different coordinate systems are given in Fig. 2, 
and have been constructed by the rotation corresponding to / and ju, of that coordinate 
system in which the function, 79 (2,y,2), has been given. The axes of this ‘‘original’”’ 
coordinate system have been chosen parallel to the corresponding axes of the coor- 
dinate system of the central site, and its origin is situated at the octahedron corner 
E on the positive z-axis of the central site coordinate system. 

Table 1 gives the coefficients for the symmetry orbitals RZ) corresponding to the 
different group elements R in the expansion’ 


Soo PE" Xoo (20) 


where P’! is the projection operator (18) and 7, is the above-mentioned symmetry 
orbital of type (19). According to (18) these coefficients are the diagonal matrix 
elements of the different irreducible representations of the inverse group element 
R“, since R is unitary. 


1 The complete basic orbital, ¢/%, is a superposition of orbitals of type (20). 


or 
bo 
— 


Pou 


R Xoo in the expression 


P19 | Tu 


H, |\B, 


Bi gt for different 7’s and i’s 


Aoy 


Agg 


(20) according to (18). 


Ayu 


ie 


é {1 |a] 1 | 23) 7 3] 12 sea a pee eee 


j 


—1 


Indices 
for the 
symmetry 
orbitals 

=R Xoo 


% 


CT ee ee ee ee a 
Say = SS a ee 
Sot et 2S Oe a ea eee 
aT AN na othind ind bs here eae wR 
TS eS SS ae ee 
gig kl SO) ONO SS Si te SP SS SS 
SN ethene cee ec 
Se Steet OO (So O'S SO OS SiO SO Sg atard 
LL ieee SS SSS See 
pe he ee ee ee 
fe Wa ie) ead SE BSI) IS) Ce) SIS) ed 
ne Bee Oo SS eee See 
a At 2 a rte eee 
ee ee ee 
Ls oe A oe oe = | 
Se eee ee 
Ss = = §— = | SS eS ee ee ee 


= = =§— §— = SS 2S 2S Se ee re Ree 


Group 
R 
of O, (7) 


Table 1. Diagonal matrix elements /(R~'), for the inverse group elements R~* of 
element 


O,,. These are also coefficients for the symmetry orbitals Z% 
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So Tee ae nee roe fe ee 
CLO 0S Se Oe ee ee 
iad acd ed ee TR ae eg em 
dee a deme me TE Seay eg hy aaa 
De I ee re ee ee a ee eee 
Sd ad elie ad ad ape) oe we ae cu apc 
SO ty ae eg es 
eh Teri ee ee eee 
Mit ake ae EM tak c ea ! 
OO le ea ee a ee ee 
ee Ee se x ce 
Se ee | 
AM SOHN ASS 
Ss se ese ewes ee re oococlcUclUlcUClULrRn hm SS 
| Tehri sa |) elie ee lee tee 
a ete ted WO ene ede i 
lo ce ce A ee | an ie 


— oe 
Wn ee og a aa Pa ae 


Ue ee ee ee ee oe | 


* ey ee s gp SEU SEN) NUCL et GUS ree 
re 8S  aQnanad Nn Worcs ri im U8 ON Gog toa Sof 8 8 = = S 8 Pee eS Se Se S&S & 
ao = © Se Sie a) Fa aa) eee i) Sa SS, See 

i ee | 
SNANANSCCAHMAMAOCMAMANOCMHMAHONANCONOCNMHAMHAO 
SON N TF MNCOCON RFI MAND ONCHEMAMAENONNANSCOMAOANAN OO & 

= = = SS = SS = = | = SS eR Le A ce ee ee ee | 

. . A STAAIndt BF MAHAN 

R SA +AANtT BA mMAFAD —§ Nad Soot ang 3 np ae a 2 2 See eee se & & 
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1 


(1) Notations, see e.g. [8] (cf. Fig. 1). 
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Table 1 (cont.) 


Group Indices 4(R~')y for different j’s and 7’s 

1 for the 

peament symmetry] R7} 
R bital Thu GF \AitAul4el4en Be [Bu | Tie | Pe | Tap | Pow Tog | Tou 

of O, (2) orbitals 

~=R Xoo Lo Se ae eee aE 1. is in AF cae £17, ih ol sneer ee Oe 

Si. 23x | Se Liat ta. 0 0108010 0. 1) 0. O-1M0 “01. 02051 
Sa, Diex Si. 1 }=—1|/-1) 1 0° 0/-0 O} 0 0-1) O0-0=1) 0 O=1).0 0 4 
3 

Say, 102 x Sty PE ee Ol eie Ol Ot Ol O— ls OO et 6 
Say 120x yi Lal jt el ete Oost OlOr laos OI Om. O 
S. LS WS Sax beatae OO ta Om OL On Omi On O 
Saz UNAS iS TeHeet Vente te | tela tet ONO Le OF Ol OF Ol LeOr 0 
6. Pblsx: C,.. Tee | le ae le ORO OF Ole) 602 0 —la0 0 
oe Tao x S,,.. Te ete Lele settee Of 1a Or OMe ORO ls Om 
o,, 100 x 9, Ey =1 j=l] 2 De ae Oe OO te DEO) 20}, 0 Leo 
oS ce ss lil a ese st ec elaine i 122 x Ovy 1 }-1}-!| 1 Act — lead OO eee O) a Le] 


(1) Notations, see e.g. [8] (cf. et Ae ht ca reer ee a) 1). 


Let us now denote by z a local operator, which means rotation of 2/2 about the 
‘“Jocal” z-axis, i.e. the z-axis belonging to the actual one of the corners Ai Bere ats)? 
the B, octahedron given in Fig. 2. Let us further introduce the local operators, 


o,=1+x, of =) =x, 


21 
0. = >. Oug 2? (a=1,2,... 6; 6 =0, 1, 2, 3). 
B 


Here, as before, x means reflexion in the yz plane belonging to the actual corner, 
and the coefficients C,, are given in Table 2. 


Table 2. The coefficients Cz, of (21). 


1 1 i! 1 1 
2 1p Wail 
3 1 0 1 0 
4 0 1 0 1 
5 1 (teak 0 
6 0 1 (1 a! 


As we see from Table 1 and the equations (18)-(20) the site functions foo are found 
as sums over atomic orbitals, localized to the different octahedron corners. In carrying 
out the summations, it is seen that all these corner localized functions can also be 
found by application of one of the operators + 0,6,, +0, _ or zero to the functions 
Yoo Xoo aNd Y1,0(4 = 0,1,2,3) respectively. In Table 3 these operators and the corre- 
sponding operands x are given for the different basic orbitals represented by the 
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Table 3. Local operators for certain local orbitals found from the starting function 
Xo: The projection of Zp) on a basic orbital, 6”, is found as the sum over all the 


transformed local orbitals for different corners A, B,... F. 
Local operators for certain orbitals x of the resp. corners 

Trreda- | pow : noe: 

cible suf- Seen eee 
sonte. | fe | Yao | Maw | Raw | aw J re 
tiong | ° TR ees ee 
Aig 0,(1+x)| 0, (1+x)} 0, (1+x)} 0, (1+x)| 0, (1+x)}] 0, (1+x) 

lu 0, (1—x) 0, (1—x) 0, (1—x) 0, (1—x) 0, (1—x) 0, (1—x) 
Aggy — 0, (1+x)}| 0, (1+x)]| —0O,(1+x)}| O,(1+x)] 0,(1+x)] 0, (1+x) 
Asy —0,(1-x)| 0, (1-—x)} —0,(1-x)} 0,(1-x)| 0,(1-x)] 0, (1-x) 
E, 1 O. (1+x)|—O,(1+x)]| 0, (1+x)|—0O,(1+x)| 03, (1+x)| 0, (1+x) 
E, 2 | —03 (1+x)] —O, (1+x)] —0, (1+x) | —0O,(1+x)| 0, (1+x)| 0, (1+x) 
EH, 1 0, (1—x)| —0,(1—x)] 0, (1-—x)| —0,(1-—x)] 0, (1-x)|] 0, (1—x) 
Ey 2 | —0, (1-x)| — 0, (1—x)| —0, (1=x) | — 0, (1—x) OS \(1— x) (Ope ee 
T's 1 0 — 0, (1+x) 0 0, (1+x)} 0, (1+x)] —O, (1+ x) 
Lae 2 0; (1—-x) 0 —0, (1—x) 0 —0,; (1-x)| -0,; (1—x) 

lo 3 0 0 0 0 0, (1—x) | —0, (1-x) 
Ty. 1 0 — 0, (1—x) 0 0,(1—x)}| 0, (1—x)| —0, (1—x) 
fhe, 2 0; (1+x) 0 —0, (1+x) 0 0;,(1+x)} 0; (1+x) 
Oho e 3 0 0 0 0 0, (1+x)| —0, (1+x) 
Lee 1 0 0, (1+x) 0 —0,(1+x)| 0, (1+x)] —90, (1+x) 

29 2 | -0, (1-x) 0 0; (1-x) 0 0; (1-x)| 0; (1—x) 

29 3 0 0 0 0 0, (1—x) | —0, (1—-x) 
Pou 1 0 0, (1—x) 0 —0,(1-x)}| 0, (1—x)| —0, (1-x) 
Lio 2 | -0, (1+x) 0 0; (1+x) 0 0, (1+x)}] 0, (1+x) 
T ou 3 0 0 0 0 0, (1+x) | —0, (1+x) 


irreducible representation and a row suffix i. Table 4 shows the effect of working 
with the different operators occurring in (21) on various functions F (x,y,z), which 
give the angular parts of the LZ and M shell orbitals. The result will be unchanged, 
if F is multiplied by a spherically symmetrical function, S(r), where r is the local 
radius vector starting from the actual octahedron corner. Introducing for the orbitals 
7% in the head of Table 3 the linear combination (19) of s, p and d atomic orbitals 
(which have the same symmetry types as the operands in Table 4), we now easily 
find the projection of 7) on the basic site orbitals. Repeating the procedure, given 
by Tables 1-3, for the symmetry orbitals localized to all the other corners of the Bg 
octahedron, we finally get the complete basic site orbitals of spd type, which are given 
in Table 5. We have previously denoted these functions by ¢7', but here they are 
characterized by the ordinary notation of the corresponding irreducible representa- 
tion, with the difference that capital letters have been replaced by lower case letters. 


524 


ARKIV FOR FYSIK. Bd 14 nr 33 


Table 4. The effect of operating with the local operators (21) on certain functions, F. 


Operator 0 | Function 0 F 
1 | 1 | a | y | z | ee | 2x | 2y | xy | LNT 
x 1 —2 Yy Zz a Sih zy —2y oe 
1+x 2 0 2y 22 2 2? 0 2zy 0 Bits aya) 
aX 0 2a 0 0 0 22% 0 2xy 0 
Z 1 y Si z 2s zy Sn VT, sei bee") 
z? 1 ae —% Zz A 09 2a xy Hp) fe 
0,=1+ 2? 2 0 0 22 aig 0 0 2xy 2 (oe — aj) 
0;,=1-2? 0 20 2y 0 0 222% Qzy 0 0 
0,=2 0, Pe 0 0 22 2:22 0 0 —2 ey —2 (a? — y?) 
0,=2 0; 0 2y —22 0 0 22Yy |= 22x 0 0 
0,=0,+90, 4 0 0 4z 42 0 0 0 0 
0,= 05-04 0 0 0 0 0 0 0 4axy 4 (x? —y?) 
0, (1+x) 8 0 0 8z 8 2? 0 0 0 0 
0; (1—x) 0 0 0 0 0 0 0 0 0 
0, (1+x) 0 0 0 0 0 0 0 8 (227 —97) 
0, (1—x) 0 0 0 0 0 0 0 Say 0 
0; (1+x) 4 0 0 4z 42° 0 0 0 4 (x2? —y?) 
0, (1—x) 0 0 0 0 0 0 0 4ay 0 
0, (1+x) 4 0 0 4z 422 0 0 0 —4 (z*—y?*) 
0, (1—x) 0 0 0 0 0 0 0 —4ay 0 
0; (1+x) 0 0 4y 0 0 0 4zy 0 0 
0; (1-x) 0 4a 0 0 0) 420 0 0 0 
O, (1+x) 0 0 Sao 0 0 0 —422 0 0 
On (Lx) 0 4y 0 0 0 4zy 0 0 0 


Furthermore, the index i = 1,2, ... 1; does not enter in this notation, since for different 
i’s the orbitals differ only by a space rotation. Thus, they have the space degeneracy 
l,. Notations such as e.g. (2s, 2p.), which occur in Table 5, signify linear combinations 
with arbitrary coefficients of the atomic orbitals inside the brackets. However, the 
coefficients of equivalent atomic orbitals localized to different octahedron corners 
A, B, ... F are common. The degrees of those parts of the approximately factorized 
secular equation (see section 3), which correspond to the basic site orbitals ¢7' in 
the L and LM approximations resp. are given to the right in Table 5. We note that 
the projections on a,,, of all the orbitals of type (19) are identically zero. Specially, 
in the L approximation also the projections of 7 on the spaces of types dgg, Ag, and e, 
are identically zero. The symmetries of the non-vanishing basic orbitals in the LM 
approximation are given by Fig. 3. In the L approximation this result is the same as 
that of Longuet-Higgins and de V. Roberts [3]. The orbitals may differ only by a 
similarity transformation. 

In the next section, we shall examine an approximate way of factorizing the secular 
equation by use of the basic orbitals. 
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Table 5. Transformed local orbitals derived from a linear combination of symmetry 
orbitals as the starting function. 


i i js 1 site are found as sums over all the 
The complete basic orbitals $?' of spd type for the centra 
seas 4 B, ... F. The dimensions of the approximately factorized secular equations are also 
given for the L and LM approximations, resp. 


Dimension of 


Basic Transformed local orbitals for the resp. corners (*) the secular 
orbital a equ. 
1 
i 4» |: Bw |e gOoe ee ee eee Lo jee 
ay S S S S S 2 5 
g 
a, 0 0 0 0 0 0 0 0 
Uu 
4 = 1 
Aq — Xx x x x AG 0 
ay = Me ay e's ry Y 14 0 1 
u 
i 2X Xo S. 2X PEK Teh Gaeisy | OC aps) 2 6 
&9 ‘ie —(X+8) | -2X | -(X+8)| -2X | X+S | X+8 2 6 
pal DAY v4 2Y Y vy ¥; 0 1 
eu | 2 ay —2Y 4 ='2:¥a| ae ¥ “6 0 1 
{ 1 0 x 0 ait y —y 1 3 
tig 2 x 0 — th 0 x x 1 3 
| 3 y 7] 7] y 0 0 1 3 
{ 1 S aay) = 7] x al 3 8 
tan 2 y S =i) sis} y 7] 3 8 
| 3 x x x x S 15) 3 8 
{ 1 x ats = HG x Yy ay 1 4 
tog 2 eae v7 x oY x x 1 4 
; | 3 y —Y 7] 4 ay =e 1 4 
{ 1 Se 7] x -—y ze ae 1 4 
ae \ 2 =9 x y —x y y 1 & 
| 3 x She x Se x PG 1 4 
Dimension of the total secular equation 24 78 


(*) Notations: S = (2s, 2p,, 3s, 3p,, 3d,.); X = (8d p._ y2)5 Y= (3d,,);% = (2p, 3p, 3d,,); y = (2p, 
3p,, 3d,,). The coefficients Cy of (19) occurring in these hybridizations are equal for equivalent 
local atomic orbitals of a certain basic orbital. The hybrids 2 and y in the table are equivalent. 


3. An approximate factorization of the secular equation 


The basic site orbitals, 6’', which have been found for the central site in the preced- 
ing section, are automatically orthogonal to each other. Since the group operators 
commute with the Hamiltonian operator of the system, the energy integrals connect- 
ing different basic orbitals of the same site vanish identically. However, the secular 
equation (15) involves integrals of orbitals localized to different sites. In the LM 
approximation the degree of this secular equation for every point k in the wave vector 
space is 78, in the L approximation this degree is 24 (cf. Table 5). Thus, in order to 
facilitate the solution of the secular equation, we shall show that matrix elements 
belonging to different basic orbitals vanish, at least approximately, also in the case 
when two different sites are involved. 
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Fig. 3. Symmetry shapes of the basic site orbitals. 


In order to do this, we consider the group operator, R, of O,, which is referred to 
the central site (see Fig. 1). We denote by R, the corresponding operator referred to 
the site n. This will be found by a translative movement of the axis of rotation (or 
the plane of reflexion or the point of inversion) of R through the distance n. Let us 
now apply this translation to an arbitrary function, (r), localized to the central 
site. Then we get the corresponding function, ¢,(r), localized to the site n, in accord- 
ance with (1). In the case of MBg, it is now readily seen that 


Rdaz = Rundm, (22) 


where m = Rn is the position vector of the site, found by the operation of R on n. 
In words, (22) means that when a group element belonging to the central site operates 
on a function ¢, localized to the site n the resulting function is also found by operating 
on the “‘translated” function ¢,, with the corresponding group element localized to 
this site. 

Since R is a unitary operator commuting with the effective Hamiltonian operator, 


h, we also have the connection 
RtAR=h. (23) 
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The basic orbitals, 42, of the different sites, n, are given by Table 5 and (1). Let us: 
now use these orbitals to construct crystal orbitals, yz, according to (3). We can then 
form the matrix elements 


[yl byl! dew deem f gilt hgh! dr. (24) 


The integrals on the right-hand side of (24) contain the positions of the sites only 
as the difference, m =n’ —n. By using this property and the relation 


we = ete, (25) 
it is easily found that 
[lt h yh! dr=du{ ve hye! dt. (26) 


Thus, the secular equation is factorized for different values of k. For equal basic 
orbitals 4”, ie. 7 =7’ and 1 =I’, (26) is also an immediate consequence of (3), since 
the corresponding Bloch orbitals form a basis for the irreducible representations of 
the group of unit translations. We have also already anticipated this result in the 
derivation of the secular equation (15). 

The elements of (26), diagonal in k, are given by a relation quite analogous to (10). 
Let us consider the energy integrals occurring in the numerator of this relation, 


ha={ op" heh" drs (27) 


Here, for the sake of simplicity, the four basic indices are omitted in the term on 
the left side. Let us introduce the symmetrical and anti-symmetrical parts of h,: 


Sa=} (hat h_a); 
(28) 
A,=}4 (ha —h_,). 
By use of (28) the numerator in the relation analogous to (10) can be written: 
> e™ha=> {Sa cos (kn)+7 A, sin (kn)}. (29) 


Now, we want to examine whether (29) is real or purely imaginary. As we shall soon 
see, no other case occurs. 


In order to make this investigation, let us consider the inversion, R =i, in O,. 
The elements of the irreducible representations for this special group element are 
1 i) =(— 1)? Ox. (30) 


Here, p is the parity of the corresponding basic orbital. The g orbitals of Table 5 and 
Fig. 3 have p= +1 and the wu orbitals have p = — 1. According to (16) and (30) 
we now find the symmetry properties 


id” =(—1)?6" (31) 


which of course are valid also if the inversion operator and the basic orbitals are 
referred to an arbitrary site, n. Thus, using (22), we have 


1 When the site index, n, is omitted, i! is referred to the central site (ef. (1)). 
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iponisel’ =(=1)? of"; (32) 


‘since i(—n) =n. On account of the unitarity property, it =i-' =i, and by use of (23), 
(31) and (32), we now have 


[o™agtidr=[(ig")*higtk dr 
=(—1)?*?'{ gh gh" dr. (33) 
Using (27) and (33), (28) can be written 
| S.=4(1+(—1)*"]ha; = An=4[1—(—1)?*?'] ha. (34) 


As is seen from (34), S, vanishes, if the participating basic orbitals have different 
parities. In the case of equal parities, A, vanishes. Since we have used only real basic 
orbitals (cf. the preceding section), h, is always real. Thus, (29) is either real or purely 
imaginary. This property simplifies the secular equation (15) considerably, since 
it might be very troublesome to solve this equation, when the elements of the secular 
determinant are arbitrary complex numbers. We now signify such an element as 
a “‘secular element’’. Let us consider three different cases: 

_ (i) The secular element contains orbitals, belonging to equal basic site orbitals. 
In this case, the parity is the same for both the orbitals, and A, vanishes. 

(ii) The secular element contains different basic orbitals, and p= p’. In this case, 
S, vanishes, and the secular element is purely imaginary. 

(iii). The secular element contains different basic orbitals, but p = p’. In this case, 
A, vanishes, but S, = h, does not generally vanish, except for the central site. Here, 
the secular element is real. 

In order to look further into the cases (ii) and (iii), we will rearrange the terms of 
the sum on the left side of (29) by use of operation on the sites n with an arbitrary 
element R~‘ of the group O, referred to the central site, according to: 


Teh = Tete gM hgsiadt 


1 


RI 


ae De he Rev | gM hes "dt. ca 
Ro kk’ 


Here the summation over R together with division by G' plays the réle of a unity 
operator, and use has been made of (23), (22) and (16). Now, if exp [ik(R ‘n)] always 
were independent of R, it would appear from (17) that (35) contains the Kronecker 
symbols 6, and 6, as factors. For the first term in the sum over the sites (n = 0), 
it is obviously true that this exponential function really is independent of R. That 
the Kronecker symbols occur as factors in this case is also a direct consequence of d’' 
forming a basis for the irreducible representations of the symmetry group O,,, referred 
to the central site. However, for the other terms it does not hold for every k. But, 
it is certainly fulfilled for all the sites, n, at the following special points in the wave 


vector space?: 
k=0; k=2(+1,+1,2+)), (36) 


1 For these k vectors, exp [i k (Rn)]= +1, independently of R, i.e. the point group O, leaves 
this exponential function invariant. 
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i.e. for the end points of the energy bands in the corner centred directions. As a matter 
of fact, when a group operator of O, is applied to n, one gets the position vector of | 
another site with the same magnitudes of its components. These are only permuted | 
between themselves. Thus, on account of the periodic properties of the exponential | 
function occurring in (35), this exponential function is independent of R at the end 
points of the energy bands in the corner centred directions of the k space. As is: 
seen, this also means that the crystal orbitals y, are symmetry adapted orbitals of | 
O, at these points, which is a well-known property. On the contrary, there is no 
other direction in the wave vector space for which the secular equation (15) is prop- | 
erly factorized by use of the basic orbitals, di. Thus, for the other directions in the 
k space the point group representations will be mixed, but probably not too much, 
In the vicinity of the points (36) the secular equation is properly factorized, yielding 
one simplified secular equation for every basic site orbital. Thus, there is rather good 
reason to believe that even the widths of the bands are only little disturbed by the 
secular elements in the cases (ii) and (iii). A further reason for believing this is that 
many of these elements are in fact still zero, and that all the non-vanishing secular 
elements in these cases are generally much smaller than in the case (i) for arbitrary 
values of k, at least in the Z approximation for MB,. Thus, we might neglect the 
secular elements in the cases (ii) and (iii). As a matter of fact, there are many other 
inevitable approximations in the present investigation, which seem to be much more 
dubious than this one. This approximation gives an immediate way of factorizing the 
secular equation (15) for every value of k. We will then get one secular equation for 
every hybrid basic site orbital given in Table 5. The problem of solving these equa- 
tions is essentially simplified on account of the fact, that all the secular elements 
in these equations are real. 

The approximate factorization of (15) can be made mainly on account of the rela- 
tion (22), which is valid for MB,, since the cubic cell belongs to the same symmetry 
group, O,, as the B, octahedron. This is not true for crystals in general. 


4. The effective one-electron Hamiltonian operator for MB, 


The experimentally determined distances in the crystals of type MB, [1] are 
rather little dependent on the different metal atoms, M, which occur. This is consistent 
with the fact that an M atom has a considerable space available in one of the “holes” 
between the B, octahedra. Hence the stability of the MB, crystal probably depends 
mainly on the boron framework, and, approximately the M atoms might be considered 
only as point charges. Furthermore, we shall suppose that the bivalent M atoms are 
doubly ionized. For the numerical treatment the case M = Ca has been chosen. Thus, 
the total Hamiltonian operator of the system under consideration can be written 


H=5(-48)-37" +45 5 ta(2e-S4) 4 


i Tip, By, B+B, 'B,B i "ip 


Wi ; (By) 1 7, _Z, - 
+33 Zu ( : - =) 48,445 spare Les (37) 
B, M Tim 


TBM i M’ M’=M’ TyweM” 


+423 


il 
+i Vij 


Here, atomic units have been used. The suffixes B and M refer to the respective 
kinds of nuclei, and the suffixes i and j to the valence electrons of the boron frame- 
work. The notation (By) above a summation sign signifies that the sum runs over all 
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the valence electrons at the boron atom B,. Certainly, the physical idea of referring 
electrons to certain atoms in the crystal cannot really be applied to molecular or 
crystal problems in general. But, for this special case it might be done mathematically 
only as a formal trick in order to divide up the total energy per electron, which is 
possible, since all the B atoms in the crystal are completely equivalent to each 
other (cf. (46)). 
4g =3 is the nuclear charge of each B nucleus unified with its two Is electrons. 
Introducing / as the formal degree of ionization, the formal charge of the M ion will 
be Zy = 2A, where 0<A<1. Complete ionization thus corresponds to 4=1. For 
the borides with three-valent M atoms, however, complete ionization would corre- 
spond to 4=1.5. The mean value of the number of valence electrons at every B 
atom is equal to ng =3+A//3. The term #, stands for the ionization energy of all 
the M atoms in the crystal. The constant potential energy of the nuclear charge points 
has also been added, in order to get the total energy of a lattice with fixed nuclei. 
The one-electron operator would correspond to the energy of one electron in 
the average field of the rest of the crystal. To obtain the total energy (37) as a sum 
over formal one-electron energies we shall also add a constant term including nuclear 
energy, but calculated per electron. This term is equal to the energy of one B nucleus 
unified with its two 1s electrons in interaction with the whole system except its own 
valence electrons, divided by its number of valence electrons, nz. We note, that it 
is not meaningless to divide up this energy, on account of the equivalency of all the 
B atoms in the crystal. The constant ionization energy, J (A), for the degree A of ioniza- 
tion of one M atom, is computed for Ca by use of experimental data [9]. This has also 
been portioned out per valence electron of boron and added to the one-electron opera- 
tor. The resulting operator is: 


(24) Zot 8) 9 yt 1 (24) VAN (8) D) \ l 
=-1} 2 + Boop + I+ er). (38 
h 2A 2 Ts 2 lyr 1+A/9 Pies 2 ca | 18 +2 ri er) ) 


Effective Z values occur also in the added constant terms of (38), since these give the 
interaction between a B nucleus and the other B nuclei together with their electron 
clouds. The sums over B and M have been performed for the 24 nearest B neighbours 
and the 8 nearest M neighbours. With increasing order of neighbours the effective 
charges tend indefinitely towards —//3 (cf. Table 6 of section 5). 

The ionic integrals connecting neighbours of the largest distance taken into ac- 
count have almost the same values as the corresponding terms calculated by using 
point charges. For the 2s orbitals these values are identical in eight figures. This is 
a reasonable argument for breaking off the sums here (at least for the L shell orbitals). 
Then, in (38) the remaining energy has been summed up as an electrostatic energy of 
point charges. In this term, eg is the potential energy of one Bg ” octahedron in the 
field of those positive and negative ions far apart, which have not been included in 
the above-mentioned sums. At very large distances the resulting charge of every B 
ion has for simplicity been considered as placed at the corresponding octahedron centre. 
Then the contribution to eg from concentric cubic shells have been computed sepa- 
rately. Certainly, this expansion might be calculated to give any value depending 
on the order of summation, since both the positive and negative terms of the expan- 
sion yield infinite non-convergent series. However, the total value of this electrostatic 
energy calculated per electron is very small. As a matter of fact the value calculated 
by use of the Madelung constant gives about — 0.03 a.u. per electron. The corre- 
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sponding value found by the present method is about —0.09 a.u. Probably the 
latter is a more realistic value, since the octahedron point charges ought to be split 
up in six different point charges, giving an increase of the magnitude of the electro- 
static energy per electron in the last term, since the electrostatic interaction within 
the site is not counted here. A comparison with the energy bands of the bonding 
site orbitals given in the tables and figures of section 8, shows that the constant 
contribution of this electrostatic energy does not change the bonding properties 
considerably. The ionization energy, J, has also the same order of magnitude. 

All the constant terms, added to the one-electron operator!, have certainly no in- 
fluence on the electron distribution, computed by use of the variation principle (14). 
The energy differences between the various states are also independent of these 
terms. However, these additional terms will give some information about the stability 
of the total crystal. 

In qualitatively comparing the eigenvalues, ¢, of (38) with those of the total operator 
(37) portioned out “‘per electron’’, €, we see that ¢ ~ €. Since € < 0 would correspond 
to a stable state at the absolute temperature zero, we may conclude that the crystal 
is stable if e <0, that is, if the eigenvalues of (38) in the ground state are negative. 
In a qualitative treatment the mean value of —e taken over all the electrons in 
the ground state might thus be taken as a rough measure of the stability of the 
crystal. : 

On account of this approximate way of choosing the zero level of the energy, only 
the difference between two eigenvalues of (38) may be given a more direct physical 
interpretation. However, the results of the present investigation show that the zero 
level of the energy is chosen remarkably well, since all the electrons really occupy 
bonding orbitals, i.e. orbitals corresponding to negative energy eigenvalues, while 
most of the unoccupied orbitals are antibonding (see Fig. 6 in section 8). 


5. Construction of the secular elements 


The basic site orbitals, 62, found in Table 5, are analogous to (11), and the corre- 
sponding secular equations of the type (15) can be solved separately according to 
section 3. As is seen in Table 5, every gz” in the expansion which corresponds to 
(11) is expressed in terms of atomic orbitals, wi, localized to the different B atoms 
of the site n. If we omit the orbital indices j, i and y, we can write 


(n) 


Pare PaB> (39) 


where the sum is taken over those B atoms which belong to the site n. The energy 
integrals occurring in (10) can then be written (with F =h): 


(m) (m+n) 
ha= > > <B’|A|B’). (40) 


Here <B’|h|B’”> is the energy integral containing the orbitals y's, and ii” 


m,B’> m+n,B”+ 


The orbital indices are omitted also here. Consequently, for convenience, we can drop 


1 In the actual computations all the added terms are not exactly constant, since the corre- 


sponding effective charges have been given a certain small amount of orbital dependence (see 
the next section). 
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the indices »’ and »”’ in the secular energy elements corresponding to (13). Thus, 


we have 


a 
n 


hy => e™ha=>. cos (Kn) ha, (41) 


since h, is real according to case (i) of section 3. 
Now, let ¢ =o * be a reflexion operator, yielding reflexion in one of the coordinate 
planes of the central site. Then, according to (22) and (23) 


=<0|h|on>=<0|oho|on>=<0|oho,|n>=hy (42) 


when both the participating orbitals belong to the same basic orbital. The last step 
follows, since the non-diagonal matrix elements of 6 vanish, and thus of=itd 
according to (16) and Table 1. The same thing is valid for the inversion and for the 
operators C,, C, and C, of O,, which signify half round rotation about the respective 
axis. All these operators have the basic site orbitals as eigenfunctions with the only 
eigenvalues + 1 (cf. Table 1). Furthermore, the basic orbitals given in Table 5 always 
have such a symmetry that at least two coordinate axes are equivalent. Taking 
these as the y and z directions and using (42), we can write the expansion (41) 


hy =heo,0,0) + 2 cos kz ha,o,o) + 2(cos ky + cos kz) he,1,0) + 
-— 2 [cos (k, + ky) + cos (k, — ky) + cos (kz +&,)-+ cos ( kz)]ha,1,0) + 
+ 2 [cos (ky + kz)+ cos (ky —kz)] hory+..., (48) 


where the components of n are explicitly written. Terms connecting sites at larger 
distances than those occurring in (43) have been neglected in the secular equations. 
However, h.o,0,0) includes some contribution also from more distant sites. The secular 
overlap elements, A, (orbital indices omitted), are completely analogous to (43). 
For orbitals of the M shell the terms consisting of overlap integrals of the highest 
order of site neighbours taken into account are very considerable, as is seen in Fig. 4. 

Introducing (38) into (43) via (40) we get a lot of integrals of different kinds. Let 
us, for brevity, introduce certain special notations for some of these integrals: 


One-centre integrals: Two-centre integrals: 
Ay =<By | B>; 
Ay =<Bo| By); By =<B,|1/rs,+1/rs| BD; 
By =<Bgll/rail Bo: C, =<B,| —4A|B); (44) 
C, =<B,| —4A| By. D, =<B,| 1/rs| By); 
Dm =<Bo|1/Tm| Bo>. J 


The index 6 (or m) for the two-centre integrals stands for the relative position vector 
between the centres of the two participating atoms. The Mulliken [10] approximation 
has been used for the three-centre integrals: 


nih - 1) <B'|B"). (45) 


BB’ TRB 


A 
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Over/aQp 
integral 


10 (Number of neighbours) 


Fig. 4. Some overlap integrals given as functions of half the distance between the two centres. 
For the nearest neighbours in the crystal the actual half distances and the corresponding number 
of neighbours have been added to the diagram. 


This approximation is rather good for small overlap integrals. Integrals of the type 
(45) occur in sums over B apart from B’ and B”. Since all the B atoms are equivalent 
in the crystal lattice, we also have the identity, 


1 1 1 ] 1 
50 (me tie =2e a (46) 


B+B,7BB, 7B’B” 


Thus, since the effective charges (53) connected with (45) are approximately inde- 
pendent of the individual atomic orbitals, these integrals would vanish when taken 
together with the corresponding contribution from the fixed charges in (38) for the 
case A =0. In other cases there is left a term with the factor —1+1 /(1 +4/9) = 
—A/(9 +A) for the sum of the integrals of the type (45). The corresponding thing is 
valid for the three-centre integrals containing 1/ry. Denoting the inter-atomic 
distances B)B and B)M with R, and R,, respectively, we find in our approximation: 
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hy =<By|h| Bo> 


(24) 1 A 
=C,— 2, by— >, 2,\D,— ——- 2) — 
ares” 2 ( : Pee) 


-213(D ae a) + : I+2* ep) A 7 
| eye aed ae 1s+2a' én) 4, (47) 
hy =<B,|h| B> 

=a Z 1 Qh 7 @) 7 

ei 7, B Aa : mt Ze 2 \ 

»— Zp tle trea (24 SB 4 ye ti+A en) Ao, (b+0). (48) 


The effective charges must have different values for the different integrals occurring 
in (47) and (48). For one-centre integrals of the type B, in (44) the notation Z, is 
used instead of Z§". For ionic integrals of the type D, in (44) and of the type (45) 
the notation Z, and for integrals of the type B, in (44) the notation Z, have been 
used for the effective charge. According to Fischer-Hjalmars [11], we can use 


Sy yd 
payyrees Melee toa, 2h <0 (n> Me) 


fotd/r) pdr 3° <1/r> 


for one-centre integrals, containing equal atomic orbitals p,, since B has the nuclear 
charge 5. In this expression, g(r) is an approximate, spherically symmetrical charge 
density function for the electron cloud at a certain atom Bo, consisting of all the 4 + 
4/3 electrons, which are not concerned by the atomic orbital g,, and r is the radius 
vector from B,. Thus, in order to determine these effective charges, it is necessary to 
compute only some simple one-centre two-electron integrals. According to Slater 
[12], it is convenient to introduce some additional averaged “exchange charge 
density”, but this has been neglected here. Thus, we are working in the Hartree 
approximation. The reliability of this method depends on the orbitals used very 
markedly. Here, the 1s orbitals have been supposed to be completely filled and all 
the other electrons—except the one under consideration—have been supposed to 
occupy 2s and 2p orbitals only. Using (49) the screening function, g(r), thus simply is 
found according to 


rg(r) =2(1 + fis)exp(— 2C1s 7) + sf (Ces 7) + pf (Cop 7); (50) 
where f(fr) =(1 +8 Cr +677? +4 0%r*) exp(—207). 
The coefficient s occurring in (50) is the number of 2s orbitals and p is the number 
of 2p orbitals in the electron cloud, and the total number is s + p = 2+A/3. The 


numerical constants £ occurring in (50) are the orbital exponents of the Slater atomic 
orbitals. For the orbitals used in the LM approximation these are given by 


Cos= 1.34, Cop = 1.205, [13]; on 
Cis=4.7, Can = Cap = 0.433, Coq = 0.333, [14]. ] 
The “diagonal” integrals in (49) apparently also depend on the different atomic 


orbitals y,. For “non-diagonal” one-centre integrals there has been used a mean value 
for the effective charges: 
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Zy(v', v"") = $[Zo (v") + Zo (v’’)]- (52) 


For the two-centre ionic integrals of the type D, in (44), Zg" is approximately 
independent of the different orbitals occurring in the integrals. Thus, it depends only 
on the distance between the different centres of the integrals. In this case the follow- 
ing approximation has been used: 


<Bo|o (B)/ra|Bo> _ _ 4, <Bolg (rs)| Bo> A eRe (53) 


<B,|1/rs|Be> »  3= 4B; byte Be ees 


Here 0(B) is the charge density of all the 5 +//3 electrons at B, rg is the distance 
between the point of integration and the nucleus B and g is the screening function 
given by (50) with s + p =3 +A/3. Since the atomic orbitals are localized to Bo, this 
approximation is rather good, at least for orbitals of the LZ shell. For large distances 
Z, tends towards —//3, as is seen from (53) and (50). In this case the contribution 
from (1s) orbitals in g is negligible. The expression (53) for Z§" has been used also 
for integrals of the type (45) and for the corresponding nuclear energy integrals. 

The values of Zg" are given in Table 6 for 2 = 0 and A = 1 and for different distances 
R,. The influence of the chosen ratio p/s on these values is rather small, since ¢,, and 
Cs, according to (51) are nearly equal to each other. In Table 6, s has been chosen 


Z,=5 


Table. 6. Effective charges Z' in a.u. (see eq. (49) and (53)). 


: Zo | : Z, 
Orbital Distance 
Py A=G pares aa Beery Avil 
28 1.5856 1.3540 | 3.2502 0.0398 — 0.2885 
2p 1.5054 1.2607 4.5964 0.0032 — 0.3297 
3s, 3p 1.0082 0.6765 6.0056 0.0002 SUES BES 
3d 1.0028 0.6700 || 7.8466 0.0000 — 0.3333 


equal to 2(2 +//3)/(3 +A/3). For integrals of the type B, in (44), Z3" =Z, has been 
chosen arithmetically intermediate to (52) and (53). Hence it is dependent on all the 
quantities ,, gy, and R, in this case. 


6. Calculation and transformation of the energy integrals 


The two-centre integrals in (44) have been solved by use of an automatic program 
constructed by the present author for the electronic computer Besk [15]. However, 
in this program the occurring orbitals are referred to parallel coordinate systems with 
common z axes. Thus, it is necessary to perform coordinate transformations for the 
orbitals of both the centres (see Fig. 5). Under such transformations the s functions 
are invariant and the p functions are transformed as vector components. The elements 
of the five-dimensional transformation matrix A for the transformation d (X,) = 
d(xg)A of the d orbitals are given in Table 7. These elements are expressed in the 
elements of the matrix B for the rotation Xp = Bxg of the local coordinate system 
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Table 7. Transformation elements A,; of the d orbitals dj= . 


a 


d; A;; for orthogonal 


3 
coordinate transformations, @, => Brynn, Xn denoting x, y, z. 
n=1 


a.\41 Dye ye dy | do | dy dg 

dy, |t(Bn- Bip + | BuBa~ ByBas| ByBo, ~ ByBaa| BaBys ~ BasBos 4)/3 (B2, — B,) 
+ Boo — Bo) 

oo By,Bys— Bs:Bay | By Bye + BysBy,| BuBsr+ ByBs:|ByBat+ BaBs| 1/3 B,,Bys 

den BuBis—By By | ByBs+ By,By, | By Bs + ByBy,| Bu Bss+BrBa| 3 By By, 

dy By By ByBe, | BisBe, t ByyBy,| BysByy+ ByBsy| BaP ss + BseBrs| /3 B,.By, 

d,s —4)/3(Bi + Bis-|—)/3 (By,Ba,+ |-V3 (ByBst |-V3(ByBs t+ | (3 B3s—)) 


— Bis — Bn) Sew apey apy + BysBss) i oB yeBze) 


at B. Here xg is a column matrix of the position vector with the point B as origin, 
and d is a row matrix with the five d orbitals as elements. By use of these transforma- 
tions we easily find the transformation of the energy integrals. For integrals consist- 
ing of d orbitals only, for example, we can introduce the five-dimensional quadratic 
matrix 


D ={ d' (xp)hd (xp-) drt. (54) 


The transformation formula is in this case 
D=A'DA”. (55) 


As has been pointed out earlier, the integrals have been calculated for the 24 
nearest B neighbours and the 8 nearest M neighbours of a given B atom. The trans- 
formation matrices B are similar for many of these integrals, but there are still many 
different types (see Fig. 5). The largest distance B-B, occurring in the computed 
integrals, is equal to the cell width, which is the same as the shortest distance M—M in 
the crystal. The largest distance B—M occurring in these integrals is somewhat larger. 
In the expansion ég, occurring in (38), however, interaction between charges at much 
larger distances has been taken into account, but this does only concern point 
charges. 

In Table 8 numerical values of some of the “diagonal” energy integrals of the 
types (47) and (48) are given. The directions B’-B”’,, symbolized by the index 6, 
are given by Fig. 5. The coordinate systems for the respective B centres are those 
given by Fig. 2, used for all the sites. This means, that “equally situated’ B atoms in 
different sites have parallel coordinate systems. The notation “diagonal” used in 
Table 8 concerns the different types of atomic orbitals. However, this notation is 
somewhat arbitrary, since the ‘‘diagonality” also depends on the way of choosing 
the coordinate systems. In any case, the given elements might be representative also 
for the order of magnitude of the many “non-diagonal” elements. 

From Table 8 it is seen that integrals connecting orbitals of the LZ shell are rather 
small at the largest distance B—B, which has been taken into consideration, but many 
of the M shell integrals are here still considerable. Thus, in order to get more reliable 
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0.735 0 —0.679 0.882 0 
=| 0.627 0.383 0.679 R= 0.251 
0.260 0.924 0.281 — 0.398 0.533 — 0.746 


Fig. 5. Some distances in the crystal, denoted by 6 = 1, 1’, 2, ... for the two-centre integrals (44), 
and the corresponding matrices for the coordinate transformations of each centre. The given 
matrices, B’ and B” can be used for transforming the occurring two-centre integrals into a form 
suitable for computation by Besk (cf. Table 7). For I and II the transformations of only the B 
coordinate systems are given. These are useful for ionic integrals of the type D,, in (44). 
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Table 8. Some “‘diagonal” elements, <p|h,|y> of (47) and (48). These are energy 
| integrals connecting different atoms (cf. Fig. 5). 


eek a ae 

28 — 0.663 — 0.342 — 0.342 — 0.149 — 0.053 — 0,012 
2p, ~ 0.125 —0.170 — 0.170 ~ 0.075 0 — 0.006 
2p, — 0.125 0.133 0.170 0.075 0 0.034 
2p, — 0.125 ~ 0.133 — 0.097 ~ 0.146 0.043 — 0.006 
38 = 0.115 — 0.430 — 0,430 — 0.386 — 0.344 — 0.293 
3D, — 0.124 — 0.392 — 0.392 — 0.340 0 — 0.228 
3p, — 0,124 0.035 0.392 0.340 0 ~ 0.007 
3p, ~ 0.124 — 0.035 0.323 0.221 0.060 — 0.228 
3 ds 0.131 0.148 —0.301 — 0.204 0.055 — 0.184 
3d,, 0.059 ~ 0.029 0.316 0.230 0 — 0.230 
3 dy 0.059 0.319 — 0.316 — 0,230 0 — 0.044 
3d, ~ 0.018 0.029 0.374 0.328 0.049 — 0.044 
poet 0.133 — 0,188 — 0.374 — 0.333 0.144 — 0.090 


results, it might seem that neighbours at even larger distances ought to have been 
taken into account. However, this would probably be rather meaningless, since, for 
example, the approximations in (45) and (53) are very rough for orbitals of the 
shell. 


7. Discussion of the methods 


As we have seen, there are many different approximations implied in the methods 
used in the present work. First of all, the chosen types of orbitals might be discussed. 

Certainly it is questionable whether it is applicable at all to use one-electron orbitals 
in crystals of the present type of structure. The answer to this question can only 
be given by the reasonableness of the found results. When working in the one-electron 
approximation it might be put in question if it is suitable to use Slater atomic orbitals 
localized to the different atoms in each site, i.e. the MO treatment. When using these 
orbitals it is maybe not convenient to break off the set of atomic orbitals already 
with the M shell. It is not evident indeed, that even more excited atomic orbitals 
need not be taken into account in order to describe the electron distribution of the 
total crystal. Thus maybe the set of atomic orbitals is no good starting point for the 
expansions of the site orbitals. However, it is difficult to find a better practical way 
of approximation. 

The cellular method would probably reduce the computations somewhat, but it 
has not been used in the present work. Since this method is not very accurate for the 
excited states, it would probably be of less accuracy than the method used for MB, 
in the present paper, when we are looking for the contribution also of M shell orbitals. 
Furthermore, the cellular method implies the use of a rather complicated unit cell 
for every single B atom in the MB, crystal. Von der Lage and Bethe have applied the 
cellular method to sodium [16], which also belongs to the cubic point symmetry 
group O,. These authors have constructed the eigenfunctions of the approximately 
spherically symmetrical one-electron Hamiltonian operator for the cubic unit cell 
in terms of so called “‘Kubic Harmonics”, i.e. linear combinations of spherical harmonics 
suitable for the cubic symmetry. The origin of the Kubie Harmonics is the centre of 
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the unit cell. In a very rough application of the cellular method to the present 
problem, these Kubic Harmonics might be chosen as a starting point in using the 
total cubic cell as unit cell, since the symmetry of the site structure might make it 
reasonable to use spherically symmetrical potentials with the origin at the centre for 
every cubic cell. However, the bonds between the different atoms within a site do not 
come to explicit expression by use of this method. In the present work the orbitals 
corresponding to these Kubic Harmonics are hybridized Slater atomic orbitals localized 
to the octahedron corners inside the cell (see Table 5). 

A strict application of the cellular method to close-packed hexagonal lattices, where 
there are two atoms inside every cell has recently been performed by Altmann [17]. 
He also showed that the symmetry properties of the basic orbitals for the total space 
group G’ of the crystal can be found for certain k-values! by using the irreducible 
representations of the point symmetry group, G. This is in agreement with the results 
of section 3 in the present paper. The whole computing process of [17] has also been 
programmed for an electronic computer. In the present problem it has only been 
practically possible to use an electronic computer for certain parts of the calculations. 

Probably the most dubious approximation used in the present work is the rough 
estimation of the electronic Coulomb energy by use of effective charges. However, 
the actual problem of computing all the two-electron Coulomb energy integrals and 
of taking these into consideration in the total energy variation problem would be 
an enormous task. In spite of using an electronic computer for part of the numerical 
calculations, it is scarcely likely that this might be finished within reasonable time. 

Furthermore, it is a rather bad approximation to cut off the sums over the sites 
as early as has been made in the present work, where atom neighbours up to the dis- 
tance between two metal atoms have been considered. On the other hand, it is cer- 
tainly no idea taking integrals for larger distances into account, not only since it 
would be a great deal of extra work, but also since the orbitals of the electrons belong- 
ing to the metal atoms then ought to be considered explicitly. For the L shell orbitals, 
the approximation of neglecting the distant sites probably gives a difference from 
the total sum over all the sites with only a few percents. For the M shell orbitals, 
however, the last term taken into account in the expansion is almost of the same 
order of magnitude as the whole sum of the preceding terms. This is very sensitive 
for the overlap integrals, since for certain points in the wave vector space the overlap 
matrices of the secular equations are not positively definite on account of the early 
cut off of the site neighbours for the M shell orbitals. This caused some trouble when 
using the electronic computer Besk, since the program for secular equations requires 
that the overlap matrices are positively definite. Thus the contribution of the excited 
orbitals to the electron distribution can only be found qualitatively. As is seen in 
the next section, however, the energy eigenvalues are generally not very much affected 
by this approximation. 

Finally, we will discuss the orbital exponents (51) occurring in the Slater atomic 
orbitals. According to more recent investigations the orbital exponents of the excited 
orbitals should be chosen larger. This might change the results somewhat. However, 
it does not affect the L approximation, and since the influence of the excited orbitals 
on account of the many approximations in the present work can at any rate only be 
described qualitatively, this may not be taken too seriously. Furthermore, since the 
valence electrons in the ground state of boron only occupy the Z shell, it still seems 


* These k-values are those for which the space groups @’ and G¥ are identical, where G* is 
the subgroup of @’ that leaves the function exp (? kn) occurring in the Bloch orbitals invariant. 
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to be well-grounded to use almost complete screening for the M shell orbitals. As a 
matter of fact, the results found in [2] yielded orbital exponents lying between 0.23 
and 0.36 for the excited 3d orbitals in the different kinds of metal borides (cf. the intro- 
duction). For CaB, the value 63, = 0.323 was found, which is very close to the value 
given in (51). 

The low values of these orbital exponents certainly depend upon the ionization of 
the metal atoms, since the valence electrons of the metal can be considered to occupy 
the excited M shell orbitals of boron. Thus, for these orbitals there are 44 K and L 
electrons almost completely screening a B nucleus, which has the nuclear charge 5. 
_ This corresponds to an effective charge somewhat larger than 2/3. Thus the orbital 
exponent for the 3d orbitals will be about 2/9 or a little more, in consistence with 
the results of [2]. As will be seen in the next section, the results of the present investi- 
gation also confirm that the valence electrons of the metal really can occupy bonding 
site orbitals mainly on account of the excited M shell orbitals. 


8. Electron distribution and energy bands 


The elements of the secular equations have been calculated for the following 
special points in the wave vector space: 


: _ (0, £2, 0), 

k, = (0, 0, 0), k,=(1 2, 0, 0), aa (0, 0, + zt), (56) 
+g eg ) 

k,= (Ex, = 7,0) k,=(0,ta,t2) ke=(ta,ta,+2). 


Cane AU ae A - 


Points in the k space giving identical energies and electron distribution have been 
put together under the same notation in (56). The plus or minus sign can be arbitrarily 
taken for every component. The directions of the y- and z-axes of the coordinate 
system in (56) are chosen to be the two equivalent directions for the basic orbital 
in question. The corresponding secular equations have been solved separately for 
all the nonvanishing site orbitals in the L approximation. The energy eigenvalues 
and the electron distribution of the bonding site orbitals in this approximation are 
given in Table 9 and Table 10 resp. Here the site orbitals are denoted by a preceding 
figure in addition giving the number of the root of the corresponding secular equation. 
The site orbitals are supposed to be doubly occupied on account of the two spin 
orientations and they are ordered with respect to decreasing bond strengths. But, 
as is seen from Table 9, the orbital 1t,,, has only one bonding direction in the table. 
Thus, it cannot be completely filled. Also some of the other bonding orbitals, namely 
le,, 1t,, and 2a,,, give only partly filled bands. The orbitals Itau, 2t,,, lt, and 2e, 
are antibonding in all directions. In roughly counting the ratio of the number of 
k points with negative energies and the total number of k points in Table 9 as a factor 
for the occupation numbers, we see that there is room only for about 18 bonding 
electrons of every site in the L approximation. However, on account of the ionization 
of the bivalent metal atoms, there are 20 valence electrons available for the bonding 
orbitals of every site. Thus the LZ approximation does not give so many strongly 
bonding orbitals as is necessary to describe all the 20 valence electrons of the site. 
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Table 9. Energy eigenvalues ¢, of the secular equations in the LI approximation. 


eS Gee 


Bact, Space | Degree €x for certain points in the k space Condes 
5; degen-| of the ARS 
site 8 
bital| CT@cY secular (0, 2,0) | (%, 7, 0) bands? 
tal oat equ. | (0,0, 0) | (2, 0,0) | (0,0, 2) | (2, 0,2) | (0, 7,7) | (x, 2, 2) 
1 aig 1 2 —0.815 | —0.805 | —0.805 | —0.788 | —0.788 | —0.753 
ltiu 3 3 0.060 2.67 12.8 1.83 — 0.706 1.99 x 
3 tiu 3 3 — 0.532 | —0.627 | —0.466 | —0.696 | —0.377 | — 0.491 
leg 2 2 — 0.553 0.694 | —0.082 | —0.486 | —0.531 | —0.423 x 
1 tag 3 1 0.068 | —0.333 | —0.334 | — 0.201 0.099 | —0.222 x 
2 aig 1 2 — 0.184 | —0.138 | —0.138 0.022 0.022 1.286 x 
1 fou 3 1 0.290 0.501 0.852 0.055 0.028 0.663 Totally 
2tiu 3 3 1.24 0.275 2.23 0.077 0.033 0.221 ae 
5 2 
1 tig 3 1 1.65 0.543 0.346 0.809 2.24 0.799 bonding 
2 eg 2 2 2.01 3.91 1.70 2.83 3.73 1.67 


1 When conduction bands occur, this has been marked with an x. 


Table 10. Kigenvectors of the secular equations for the bonding orbitals in the DZ 


approximation. 
PR MOAtorah Coefficients C, of the atomic orbitals for different points in the k space 
ecbilel viiay Nab: eee 
y (0, 0, 0) (zz, 0, 0) (0, 0, 2) (7030 '3t) (0, 2, 7) (7, 7, 7) 
“Be {2e 0.996 0.991 0.991 0.989 0.989 0.993 
1g |2po 0.092 0.130 0.130 0.147 0.147 0.119 
28 0.925 0.817 0.881 0.805 0.999 0.719 
1lt,,, 2po = 0-101 — 0.187 — 0.196 — 0.461 0.001 0.579 
2px — 0.367 0.546 — 0.432 — 0.374 — 0.038 — 0.384 
28 0.980 0.939 0.974 0.999 — 0.454 0.865 
3th, 2po 0.035 0.228 0.005 0.034 — 0.014 — 0.003 
2px 0.195 0.259 0.225 0.018 0.891 0.502 
re 28 0.984 0.991 0.964 0.952 0.977 0.964 
9g 2po 0.180 0.133 0.266 0.307 0.214 0.267 
lig, 2pm 1 1 1 1 1 1 
Dy | 2s 0.152 0.336 0.336 0.549 0.549 0.749 
1g 2po 0.988 0.942 0.942 0.836 0.836 0.662 


ee ee eee ee eS Sa 


(*) “Atomic orbital” here means the basic site orbitals of Table 5 with only one type of atomic 
orbital. The notation o in this table corresponds to zin Table 5 and 2 corresponds to # and y. 


In the LM approximation, however, there is a sufficient number of strongly bond- 
ing site orbitals. This is seen from Table 11 and Fig. 6, which give energy eigenvalues 
and energy bands for some bonding orbitals resp. The corresponding electron distri- 
bution is given by Table 12. The coefficients C, in Tables 10 and 12 are those of (19). 
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Table 11, Energy eigenvalues e, of the secular equations in the LM approximation. 


a a ee 
Space | Degree éx for certain points in the k space 


a degen-| of the 3 cae 
; eracy | secular (0, 7,0) | (2. 2, 0 ees 
orbital) 1, | equ. | (00,0) | (2,0, 0) | (0,0, ) | (ae, 0, x) | (0,7, 2) | (x,2,2) | bands" 


Sp ear eee se a ee ee ee ee eee 


l aig 1 5 mend (AON 13} — 0.861 — 0.861 
ltin 3 8 0.657 — 0.728 x 
Ley 2 6 — 0.57 ~ 0.38 — 0.45 (a) 
3hu 3 8 — 0.504 — 0.345 
6 tu 3 8 — 0.391 — 0.421 
3 dig il 5 — 0.382 |-— 0.388 | —0.388 
3 tog 3 3 — 0.364 0.561 =(0:273 x 
l tag 3 + 0.029 | —0.350 | —0.349 AV hel ue 
5 eg 2 6 — 0.34 — 0.284 
4g 2 6 — 0.34 =\():2 
4 tiy 3 8 —(OoL2 — 0.247 
fist 3 4 1.28 ~ 0.320 0.5 0.4 Ps 
2 tog 3 4 — 0.253 0.037 0.066 airs de oe 
Tt | 3 3 | +023 ~ 0.253 
5tiu 3 8 1.99 — 0.294 x 
6 eg 2 6 — 0.206 — 0.274 
4ay| 1 me |e O68 
8 tiu a4 8 = 0.227 0.446 x 
3 eg 2 6 ~0.179 
Ray | 1 5 | —0.161 | —0.064 | —0.064 («) 
5-3 ne} 3 ~ 0.147 ~0.160 
5 aig 1 5 0.038 | -— 0.142 | —0.142 r 
ee 4 — 0.127 
Ttiu 3 8 0.051 0.020 
Fee 3 4 0.082 
dou 1 1 5.2 0.3 0.1 
oy a 3 1.2 0.4 0.7 | Totally 
eu 2 1 0.5 anti- 
4 tou 3 4 2.9 U6 bonding 
2eg 2 6 4.2 2.5 2.9 
Ste | 3 8 3.5 2.8 
ae lie T 1 11.6 
2 tou 55) 4 13.5 


1 See the note under Table 9. The notation (2) means that the conduction bands can be seen 
only from Table 9. 


The explicit expression of the site orbitals are finally found by using Table 5. As we 
see from Table 12, the distribution of d orbital electrons is remarkably essential also 
for the most stable configurations in good agreement with the supposition in [2]. 

On account of the uncertainties discussed in the preceding section for the M shell 
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€ (a.u.) 


Ik 


Fig. 6. Regions of the energy bands for some of the bonding site orbitals. Bands for orbitals with 
large contribution from the M shells (e.g. 3ty,, 3a,, and 6¢,,,) are rather uncertain, 


orbitals, it is not possible to find the solutions of the secular equations for some 
of the directions in the k space without working out more terms in the expansion 
over the site neighbours. Hence there is some open space in the tables. However, in 
spite of the somewhat incomplete and inaccurate values for the excited atomic 
orbitals, Table 11 constitutes a good complement t9 Table 9. As is seen, the inclusion 
of excited orbitals does not disorder the energies of the orbitals in Table 9. This only 
gives a little lowering of the lowest lying levels. The main point here is, however, that 
some of the strongly bonding site orbitals, which do not excist in the L approximation, 
namely for example 6t,,, and 3a,, are stronger bonding than several of the weakest 
bonding orbitals in the Z approximation. This means, that in the LM approximation 


there is a sufficient number of strongly bonding orbitals for the 20 valence electrons 
of every site. 
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Table 12. Kigenvectors of the secular equations for the strongest bonding orbitals 
in the LM approximation. 


Coefficients C, of the atomic orbitals for different points in the k space 


Basic Atomic 
site pepiial (0, , 0 
orbital te ( ) (0, 0, 0) (x, 0, 0) (0, 0, 2) (0, It, 7) | (a, 2, 2) 
2s 0.917 0.923 0.923 
2pe 0.145 0.097 0.097 
la,, 38 — 0.173 ~ 0.107 — 0.107 
3 po 0.032 0.246 0.246 
3do ~ 0.327 ~ 0.258 — 0.258 
28 0.882 0.953 
2pa ~0.161 —0.031 
2px — 0.313 ~ 0.007 
e 38 — 0.257 — 0.134 
lu 3pa 0.041 0.185 
3 pa 0.043 — 0.000 
3do 0.161 0.125 
3d2 0.051 0.153 
28 0.948 0.796 0.916 
2po 0.282 0.073 0.325 
38 — 0.146 ~ 0,142 — 0.216 
le, 3pa — 0.005 0.505 ~ 0.044 
3da — 0.020 — 0.274 0.021 
340 — 0.013 0.106 0.074 
28 0.919 0.132 
2Qpo 0.045 — 0.255 
2px 0.321 0.584 
38 0.036 — 0.506 
3 thy 3 po 0.036 — 0.104 
3pm 0.134 —0.310 
3do — 0.107 0.409 
3dx — 0.132 0.215 
2s ~ 0.196 ~ 0.418 
2pe —0.081 0.079 
2px ~ 0.519 0.757 
38 — 0.436 ~ 0.354 
6th, 3po 0.624 0.068 
3 p27 0.252 0.293 
3do ~ 0.180 0.047 
3dz ~ 0.105 — 0.167 
28 — 0.455 — 0.337 — 0.337 
2pa 0.238 0.308 0.308 
3a, 38 0.718 0.797 0.797 
5 3po — 0.465 — 0.376 — 0.376 
3do 0.064 0.119 0.119 


() The symbols o, 7, 6 of the d orbitals correspond to the orbitals of the types 3d,:, 3d,, 
and 3d,,, 3d,+_,*, resp., in Table 5. See also the note under Table 10. 
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Furthermore, it is seen from Table 9 and Fig. 6, that the doubly degenerate orbital 
le, is strongly antibonding for the k vector (z, 0, 0), which is perpendicular to the 
plane containing the 2s and 2p orbitals (see Fig. 3). These energy bands are thus 
only partly filled, implying metallic conductivity. Certainly, these conduction bands 
will be a little disturbed on account of the fact that the secular equation is not exactly 
factorizable in this direction of the wave vector space (see section 3). But, the occur- 
rence of these conduction bands depends mainly on the fact that we are working outside 
the tight-binding approximation. Furthermore, the most distant neighbours taken 
into account only modify the energies with about one per cent. Thus, the existence of 
these bands seems to be quite certain. An interesting point in this connection is, 
that for sodium the electrons also are essentially free in the corresponding direction 
[16]. As is seen from Table 11, similarly also 1t,,, and the orbitals of type t,, are partly 
antibonding. 


9. Discussion of the results 


The results given in the previous section are partly in contradiction to the results 
found by Longuet-Higgins and de V. Roberts [3], who have found a closed shell 
structure for MB, without any conduction bands. The experimentally found metallic 
conductivities of the metal borides with bivalent metal atoms have been explained 
by these authors as effects caused by impurities or lattice defects of the crystals. 
Of course, this cannot really be disproved, but the results of the present investigation 
show a more natural way of explaining the occurrence of these conductivities, since 
several of the most strongly bonding site orbitals show up conduction bands. 

The sequence of the strongly bonding site orbitals found by the present method 
in the L approximation is mainly in accordance with the results found by [3]. In [3] 
the stability has been explained by use of only L shell orbitals. The present investiga- 
tion shows, however, that also M shell orbitals are necessary for explaining the 
stability of the crystal, since some of the site orbitals, which only exist in the LM 
approximation, are more strongly bonding than the weakest bonding orbitals in the 
LI approximation. 

The high stability of the crystal is very well explained by the present results, 
since all the occupied orbitals are strongly bonding. Also for the borides with three- 
valent metal atoms there is a sufficiently large bonding energy for the extra electron 
to compensate the ionization energy and the Coulomb repulsion energy, which come 
in addition. Thus also crystals of this type are very stable. Since in this case there 
is an odd number of valence electrons in every site, these crystals have large metallic 
conductivities in good agreement with the experimental facts. 

Table 12 shows that the M shell orbitals are very essential for many of the strongly 
bonding site orbitals. This result confirms the necessity of taking these orbitals into 
account. When the occupation numbers of the M shell orbitals are large, however, 
they can be correct only for the order of magnitude, since many of the integrals 
connecting these orbitals have barely reached their maxima for the largest distance 
oecurring in the expansion over the atomic neighbours (see Fig. 4). The presence of 
some large occupation numbers for the excited orbitals shows that it is rather doubtful 
to work with only non-excited atomic orbitals in an attempt to describe the electron 
distribution in a crystal. Also when introducing some excited orbitals, it is rather 
dubious if the present method is always acceptable, since it is too laborious to work 
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out and transform all the integrals for a sufficiently large number of neighbours, 
and also since it is rather impossible to know, if even more excited orbitals should have 
been taken into account. 

Taking into consideration all the rough methods and all the other inevitable 
approximations used in the present work it is very remarkable that the results really 
give so extremely good explanations for the properties of the crystal. 


APPENDIX 


Derivation of an expression for the effective Hamiltonian operator 


Let us consider the two-centre integrals occurring in (6): 


Jn=><kK'|H®|kk); 9K, => kk’ |H®|k’ kb). (57) 
k’ k’ 


For basic site orbitals with equal parities, the factor which corresponds to the normal- 
izing factor of (3) is real. Thus, in introducing the site orbitals according to (3) 
we can write 

e k’(n’’—n’) 


© We = > se be hs 5 
vie = > pygemy bebe (58) 


n’n” 


This expression occurs in the integrands of both J, and K,, summed over k’. Thus, 
let us introduce the notation 


V> oo k’ 


1 ek (n’’—n’) yg Sate ei ka’—n) 
Dyn = lim 25 R= A a | | | seach dk, dks. (59) 
m re 000 ss 


Especially for orthogonal site orbitals, it must be realized that 
Daca ae Ow n”’ (60) 


is a Kronecker symbol. This is also merely a consequence of the fact that the elements 
of the one-electron density matrix are invariant for a transformation (3)=(4) in the 
orthogonal case. 

By use of the overlap matrix 


A = [ & (r) @ (r) dr, (61) 


where &(r) is a row matrix with the site orbitals 4, (r) as elements, we may construct 
orthogonal site orbitals y(t) = %(t — 0) according to the Léwdin [4] formula, 
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where x is a row matrix with y, as elements. On account of (62) we see that 
Dix Xie (63) 
Here is a row matrix of the crystal orbitals, y.. Summing over k’ in (58) yields 
Yyy=O@De (64) 


where the elements of D are given by (59). Comparing (63) and (64), which can be 
considered as the unity operator in Hilbert space, we see that D =A’. Now, from 
(64) 


> Wer (r’) Wx’ (r’’) ae 2D Dara pa’ (r’) pn” (r.) = (65) 
Formally, this is an element of the charge density matrix, although the sum here 
goes over all k’. Using (65) and (3), (57) gives 
Ju=( Vp An) se Cae OS Dorag ie WW |) Aten 


‘,m’" 
” 
2 


fies (66) 
K,=( ve ae Ne} -1 8 ik (m’ Se ee <m’n’ | H™? | n’’ m’”). 


For orthogonal site orbitals, these formulae are in accordance with those given by 
Wohlfarth [18]. However, since V is the infinite number of sites in the crystal, it is 
convenient to reduce these expressions furthermore. As a matter of fact, the integrals 
in the site orbital space do not change their values, if all the centres, which occur in 
the integrals are changed by a common translation belonging to the total symmetry 
group of the crystal. Thus, it is allowed to subtract the vector m’ from all the four site 
indices of every term in (66) without changing the result. Hence, if we make the sub- 
stitutions m’’=m +m’, n’=p’+m’ and n’ =p’ +m’ into (66) we see that the 
expressions inside the sums are independent of m’. Thus, the summation over 
m’ gives a factor V, which is shortened by the factor 1/V. Hence we can write 


> em ded pS ert en 
Ix = Ser Ae ; Ly yam A? (67) 
where af ee Dy—» <0 p’ | H*?| m p’’» (68) 
and Ku= 2 Dy» <Op'|H® |p" m). (69) 


In the special case of orthogonal site orbitals there is found 
=2 (On| H”|mn> (70) 
and Ka=><On|H™"|nm). (71) 


In the general case we can introduce the following notation (cf. (6)): 
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Fi=Hi+4, Jd, —A, K,. (72) 


This gives the expression (10) of ¢,. By using (11) for the site orbitals, the equation 
(12) is formally unchanged, if we use the following expression for Fy””: 


See ee 
ote ze {(2 a eee Ont 2,, Dare 
IgA a ve ON ee eee 
Bian |e Gaal Aa atl Opie iF 


The site indices have been put below the orbital indices in the one- and two-electron 
integrals above. Since the coefficients C, occur in (73), Fy” are not independent of 
the variation in (14). 

In the Hartree approximation the Fock operator working in the gy space now can 


be written 


OKT ee 
F=H'+ Ay DC Ow & Darr | git” (2) — oh (2) dt,. (74) 


"12 


By use of (74) we find the elements, which in the one-electron approximation are 
replaced by the energy elements hy” of (15): 


Pe = Seo" a) Fo ary. (75) 


The effective one-electron Hamiltonian operator (38) corresponds to the Fock operator 
given above with suitably chosen coefficients C, and the inverse overlap matrix D 
replaced by the unity matrix. 
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